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ABSTRACT 


The  unsteady  expansion  of  gas  clouds,  bounded  by  a  vacuum,  is 
considered  for  both  continuum  and  free  molecular  flows.  Known  analytical 
solutions  are  reviewed  with  emphasis  on  the  asymptotic  nature  of  the  flow 
after  long  times.  An  approximate  analytical  solution  is  presented  for  the 
continuum  expansion,  after  long  time,  of  an  initially  uniform  gas  cloud.  For 
a  plane  expansion,  the  approximate  solution  is  in  exact  agreement  with  the 
known  exact  solution.  For  cylindrical  and  spherical  expansions,  the  ap¬ 
proximate  solution  gives  results  for  density  (at  the  axis)  which  are  about 
10  percent  higher  than  those  obtained  by  a  numerical  integration  of  the  equa¬ 
tions  of  motion.  (Only  two  numerical  solutions  are  available  for  the  latter 
comparison.)  The  continuum  and  free  molecular  expansions  of  an  initially 
uniform  gas  are  found  to  be  quite  similar,  after  long  time,  except  for  the 
flow  in  the  vicinity  of  the  leading  edge  of  the  expansion.  For  a  given  gas, 
the  density  at  the  axis  decays  more  rapidly,  with  time,  for  the  continuum 
case.  The  extension  of  these  results,  to  hypersonic  jets,  bounded  by  a  vac¬ 
uum.  is  also  noted. 
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NOMENCLATURE 


sound  speed 

hypergeometr ic  function 
y/2a^ 

Mach  number 
elemental  mass 
(y+  1 ) / 2(y -  1 ) 
pressure 

ordinate  of  leading  edge  of  expansion 
spatial  coordinate,  transverse  for  jet  flows 
time 

axial  velocity  for  jet  flow 
velocity,  transverse  for  jet  flows 
axial  coordinate  for  jet  flows 
characteristic  Riemann  invarients 
gamma  fraction 
specific  heat  ratio 
characteristic  streamline  slope 
similarity  variable 
density 

dimensional  index  (0,  1,  2) 


INTRODUCTION' 


Spate  flight  has  in’ensified  interest  in  the  unsteady  expansion  of  gas 
<  loads  into  a  vacuum  A  related  problem,  the  lateral  expansion  of  a  hyper¬ 
sonic  iet  bounded  by  a  vacuum,  is  also  of  current  interest 

Several  analytic  solutions  exist  for  unsteady  expansions  into  a  vacuum. 
Reform  es  1  and  2  have  treated  plane,  cylindrical  and  spherical  self  similar 
Hows  These  are  flows  where  time  c  an  Oe  eliminated  from  the  equations  of 
motior  reduc  ing  them  to  ordinary  differential  equations  which  can  be  inte¬ 
grated  -n  closed  form.  Only  special  initial  conditions,  for  a  gas  cloud,  can 
lead  to  self  similar  solutions  In  particular,  die  initial  density  distribution 
<a:  no*  be  uniform  The  plane  problem  of  an  initially  uniform  gas  cloud, 
expanding  irto  a  vacuum  has  been  solved  analytically  in  Reference  3  The 
solution  ior  an  initi  ally  uniform  cylindrical  or  sphen.al  gas  cloud  requires 
numerical  ir.tegiation  of  the  equations  of  motion.  The  above  studies  use  con¬ 
tinuum  equations  of  motion  If  the  mean  free  path  ir.  a  gas  cloud  is  of  the 
order  lor  larger!  of  the  cloud  diameter  the  resulting  expansion  can  be  treated 
as  a  free  molecule  flow  Reference  4  has  obtained  solutions,  in  quadrature 
form,  tor  the  tree  molecular  expansion  of  arbitrary  gas  clouds 

!n  the  present  study,  the  asymptotic  nature  of  the  flow  field,  after  long 
times,  if  examined  for  each  of  the  analytical  solutions  mentioned  above 
These  are  used  *  a  guide  to  develop  an  approximate  asymptotic  solution  for 
the  expansion  of  an  initially  uniform  cylindrical  or  spherical  gas  cloud  Ref¬ 
erence  5  has  discussed  tne  asymptotic;  natuie  of  the-  flow,  near  the  origin,  for 
the  plane  problem  of  an  initially  uniform  gas  cloud  expanding  into  a  vacuum 
The  present  study  is  a  aeralization  ot  Reference  5  in  that  the  entire  flow 
field  between  the  or;0.n  and  the  leading  edge  of  the  expansion  is  considered. 

Ir.  addition,  cylindrical  and  spherical  expansions  are  considered  The  ap- 
plic  ation  ot  the  unsteady  plane  and  cylindrical  flow  results  for  finding  the  lat¬ 
eral  spreadinc  of  a  steady  hypersonic  jet.  bounded  by  a  vacuum,  is  also 
discussed 
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UNSTEADY  EXPANSION  OF  GAS  INTO  VACUUM 


General  Considerations 

The  equations  of  motion  lor  unsteady  flows  in  one  spatial  variable,  7, 
may  be  written 

Continuity:  +  cr  =0  (la) 

9t  9r  I' 


Momentum 


Dv  _  9v  1  9p 

dT  *  X  3r  p  9r 


Yl  P 


11 

dr 


(lb) 


when  the  superscript  bar  indicates  a  dimensional  quantity  and  cr  =  0,  1,2  for 
plane,  cylindrical  and  spherical  flows,  respectively.  See  Nomenclature  for 
definition  of  the  other  variables  Eq .  lb  and  all  subsequent  developments. 

•  -V 

assume  that  the  entropy  of  the  gas  is  uniform  (i.  e.  the  relation  p  =  cp  , 
where  c  is  a  constant,  has  been  used  to  eliminate  the  variable  p  in  Eq.  lb  so 
that  there  are  now  two  dependent  variables,  p  and  v  ). 

Eq .  I  can  be  nondi  mensional  i  zed  by  introducing  the  following  variables 


r 


(2) 


in  all  subsequent  developments.  R  ^  will  be  taken  to  be  the  initial  (  t  =  0)  gas 
cloud  radius  For  initially  uniform  gases,  and  are  the  initial  density 


I 


and  sound  speed,  respectively.  For  initially  non-uniform  gases,  pQ  and 
are  evaluated  at  r  =  0 ,  t  =  0.  Substitution  of  Eqs.  Z  into  Eqs.  1  yields 


9p  v 

“SF 


+  cr 


pv 

r 


=  0 


(3a) 


9v  9v  y-Z  9p 
9T+  v  37  =  '  p  37 


(3b) 


which  are  the  non-dimensional  equations  considered  herein. 

As  pointed  out  by  Reference  6,  the  pressure  gradient  term  in  the  mo¬ 
mentum  equation  (right  side  of  Eq.  3b)  becomes  negligible,  after  long  times, 
when  studying  the  expansion  of  gases  into  a  vacuum.  Eq.  3b  can  then  be  in¬ 
tegrated  to  yield 

v  =  r /t  (4) 


Reference  6  refers  to  such  flows  as  "inertia  dominated".  That  is,  each  gas 
particle  moves  with  a  constant  velocity  and  is  unaffected  by  the  other  particles. 
Substitution  of  Eq.  4  into  Eq.  3a  shows  that  the  density  distribution  must  have 
the  form 


pt'r+l  =  f(r/t)  (5) 

where  f(r/t)  is  a  function  of  r/t.  The  precise-form  for  this  function  depends 
on  the  earlier  motion  of  the  gas  cloud.  The  motion  is  self  similar,  after  long 
times,  since  v  and  pt°"^  are  functions  only  of  r/t.  All  solutions  discussed 
herein  will  have  this  asymptotic  behavior. 

Sell  Similar  Flows 

References  1  and  Z  have  investigated  expansions  into  a  vacuum  which  are 
self  similar  at  all  times.  These  flows  will  now  be  discussed. 
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Let  R(t)  be  the  leading  edge  of  a  gas  cloud  which  is  expanding  into  a 
vacuum.  The  gas  particles  are  then  contained  in  the  region  0  —  r—  R(t).  A 
vacuum  exists  for  r>R(t).  Define  a  similarity  Variable 

p  -  r/R(t)  (6) 

so  that  p  =  1  corresponds  to  the  leading  edge  of  the  gas  cloud  and  the  gas  is 
confined  to  the  region  0—  p  —  1.  Self  similar  solutions  are  found  by  assuming 
the  dependent  variables  have  the  form  v  =  g  (t)p  (p)  and  p  =h(t)f(p),  substi¬ 
tuting  into  Eqs.  3,  and  finding  the  form  of  g  (t)  and  h  (t)  such  that  the  result¬ 
ing  equations  are  independent  of  t.  The  resulting  solution,  for  a  gas  expanding 
into  a  vacuum,  is  (Refs.  1  and  2) 


P 


=  R'(cr+1)  [  1-p2 


1 

Y-l 


(7a) 


V 


,dR. 

(_dT‘)  11 


(7b) 


whe  re 


dt  vTTT  y-1  1 


This  solution  satisfies  the  boundary  conditions: 


at 


t  ^  0;  R  r  1  , 


(8) 


(9a) 


t  ^  0, 


r  -=  o; 


P  =  1 


(9b) 


t  >  a,, 

n  =  i; 

p  =  0 

(9c) 

o' 

A  I 

ii 

O 

v  =  0 

(9d) 

v  =  0 ,  and 

p  varies 

d 

3 

at  r  =  0  to  0  at  r  =  R. 

At  later  times,  v  is  not  zero,  and  varies  linearly  with  r. 


v/(dR/dt) 

nCT+1 


The  variation  of  v/(dR/dt)  and 
pR0^  with  r)  is  illustrated  in 
Sketch  1. 


v/(dR/dt),  pR0^ 

Sketch  1 

Variation  of  v(dR/dt)  and  pRa+*  with  r| 

The  location  of  the  leading  edge  of  the  expansion,  at  any  time,  can  be 
found  by  integrating  Eq.  8,  which  can  be  put  in  the  form 


1  2 


<r+  1  v  "  1 


'1  VI  -R  "(<T+  1  >  (^_1) 


(10a) 


r  -1  d1/n 

/  sec  R  o 
Nq  /  (sec  ©)  1  +  N°  d  0 


(10b) 
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where  N0  =  2  / (tr+  1 )  ( y - 1 )  •  Eq.  10b  is  obtained  from  Eq .  10a  by  letting 
R  =  sec^o  0.  Eq.  10b  can  be  integrated  in  closed  form  when  Nq  is  an  integer. 
For  example 


-==L=  -i,.  t  r'N/R^-l  for  N  =1  (11a) 

-s/T+T  V'1  ° 

=  ^R(R-1)  +  lnf^R-1  +^R)  for  N  =  2  (lib) 

=  (RZ/3  +  2)  ^/rZ  /  3  -  1  for  N  =  3  (11c) 

'  o 

Eq.  11a  includes  the  case  a  =  2,  y  r  5/3  Eq.  lib  includes  the  cases  cr  =  1 , 
y  =  3/2  and  a  =  2 ,  y  -  4/3.  Eq.  11c  includes  the  cases  or  =  0 ,  y  =  5/3; 
cr  =  1 .  y  =  4/3  and  c r  -  2.  y  -  11/9.  Thus  Eq.  10b  permits  the  complete 
closed  form  solution  for  a  large  number  of  rases  of  interest.  The  reduction 
of  Eq.  10a  to  Eq.  10b  does  not  appear  to  have  been  presented  previously.  The 
variation  of  R  with  t  is  illustrated  in  Fig.  1  which  presents  results  for 
y  ”  4/3,  5/3  and  <y  - ■  0,  1  and  2. 

The  asymptotic  nature  of  the  flow,  after  long  times,  will  now  be  noted. 
Eq.  10a  shows  that  for  R  ^°"+  ^  ^  ^  <<  1  (i.  e.  after  long  times). 

R  =  _L - £_  t  (12) 

qfo+T  y  - 1 


Thus,  after  long  times,  the  leading  edge  moves  with  a  uniform  velocity 
dR  2 

"cTT  ~ +  V ( - jy  Eor  o'  =  0,  this  velocity  is  the  same  as  that  associated  with 


Figure  1.  Lending  Kdge  if  Fxpansion  for  Self  Similar  Flows 

( I-  rom  Ref.  1  ,  pg.  27b) 
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the*  expansion  of  an  initially  uniform  gas  cloud  :nto  a  vacuum'  However, 
for  <r  1  .  d .  the  final  velocity  is  less.  by  the  factor  i/^oTT.  than  the 
corresponding  voloc  it y  lor  an  nut; ally  uniform  gas  c  loud  The  reason  for  this 
is  that  ;n  the  present  use  the  particles  at  the  leading  edge  of  the  expansion 
have  no  energy  at  t  0  and  are  subsequently  accelerated  by  the  particles  in 
the  interior  of  the  gas  cloud  For  u  -  0.  there  are  sufficient  particles  at  and 
near  the  energy  corresponding  to  conditions  at  t  r  0.  to  drive  the  lront  to 
the  limiting  velocity  dR/dt  -  2/ty-l)  With  increase  in  <j .  there  is  a  decrease 
m  trie  relative  number  of  particles  which  have  energy  at  or  near  the  value  at 
t  r  0  Hence  the  front  is  accelerated  to  lower  velocity  as  a  increases 

Alter  long  times  the  velocity  and  density  are  uiver  by  (from  Eqs  7  and 
li) 


v 


r  <  t 


(13a! 


.  o"4-  i 


<\  -  1>^  ’ 


(1  3b) 


where  rj  y,-.  -  1  ±—~  These  expressions  arc  in  agreement  with  Eqs  -1 
and  - 

Analytical  Solution  for  r  -  0  and  Uniform  Initial  Density 

Red  3  lias  given  the  characteristics  solution  tor  the  plane  (a  -  0)  ex¬ 
pansion  into  a,  vacuum  of  a  gas  whi.ii  is  initially  ur.ilorm  This  solution  has 


'  Vor  an  initially  uniform  gas  cloud.  dR/dt  -  2/(y-l).  for  all  t  >  0.  re¬ 
gardless  of  the  value  of  <r  (The  case  ot  an  initially  uniform  gas  and  ,  :  0  is 
treated  m  the  next  sec  tion)  This  is  true  because  at  t  -  0  the  flow  in  the  vi  - 
i  mity  of  r  =  R  behaves  as  through  a  0  The  leading  edge  achieves  the 
limiting  velocity  dR/dt  --  Z/(y-l)  and  is  unaffected  by  the  subsequent  expan¬ 
sion  oi  the  remainder  of  the  gas  See  Ref  7  for  further  discussion 
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been  discussed  in  greater  detail  in  Refs.  5  and  6.  The  uniform  gas  originally 
extends  from  r  =  0  to  r  =  1  where  r  =  0  can  correspond  to  a  wall  or  to  a 
plane  of  symmetry.  The  expansion  is  assumed  to  start  at  t  =  0.  At  first  a 
simple  expansion  wave  is  generated  which  is  centered  at  r  =  1  ,  t  =  0.  This 
region  is  bounded  by  the  lines  (3  =  -1,  (3  =  1  and  a  =  1  in  Fig  2.  The  reflec¬ 
tion  of  this  wave,  at  r  =  0,  results  in  an  interaction  region,  also  indicated  in 
Fig.  2.  The  solution  for  both  regions  is  discussed  herein.  In  this  discussion 
it  is  convenient  to  use  the  local  speed  of  sound  a  =  a/a^,  instead  of  p,  as  a 
dependent  variable.  These  are  related  by,  for  isentropic  flow, 


2  /  ( y  -  1 ) 

P  =  a 


In  the  simple  wave  region,  the  dependent  variables  depend  only  on 
(r  -  l)/t.  The  solution  of  Eqs.  3  for  v  and  a  can  then  be  shown  to  be 


(15a) 


a  -  — —  f  1  -  111  I_l| 
"  v+1  1  2  t  1 


(15b) 


The  leading  edge  of  the  simple  wave  corresponds  to  a  =  0  and  its  location 
is  given  by  R  =  r  =  1  +  2t/(y-I)  .  The  leading  edge  moves  with  velocity 
dR/dt  =  2/(y-l)  as  previously  noted.  The  trailing  edge  of  the  simple  wave 
is  defined  by  the  line  r  =  (1  -  t).  The  fluid  between  this  line  and  r  -  0  is 
uniform  at  the  initial  values,  v  0 .  p  a  1 

The  solution  for  the  interaction  region  requires  the  introduction  of  the 
Riemann  variables 


y  - 1 

a  =  a  +  ~ —  v 


(16a) 


|3  -  a  -  — ,-v 


(16b) 


9 


0 


\ 


If! 


0 


The  quantities  a  and  (3  are  constant  along  characteristic  lines  which  propagate 
with  velocity  v  +  a  and  v  -  a,  respectively.  These  characteristic  lines  are 
illustrated  in  Fig.  2.  (In  the  simple  wave  region  a  =  1  while  P  varies  from 
-1  at  the  leading  edge  to  +1  ait  the  trailing  edge).  From  Eqs.  16  it  follows 
that 


v  =  7TT  <Q  -  P> 


(17a) 


a 


=  i  (u  t  p) 


(17b) 


With  these  relations,  Eqs.  3  can  be  inverted  so  that  (r,  t)  are  considered 
functions  of  (n,p).  The  resulting  Eqs.  are 


dr  1 

^  y  -  1 


n  + 


(18a) 


3r  1 


Hi 

2 


(18b) 


Eliminating  r  yields  a  single  second  order  partial  differential  equation 


a2t  n  .at  at 

+  3p 


)=  0 


(19a) 


where 


The  boundary  between  me  simple  wave  region  and  the  interaction  region 
(line  a  r  ]  in  Fig.  2)  is  defined  by  (see  Refs  3.  5  or  6) 

MI.PIMt^)N  WO) 

where  0  varies  from  I  to  - 1  as  t  goes  from  I  to  on  An  explicit  expression 
for  r.(a .  in  the  interaction  region  is  found  by  integrating  Eq.  19a,  with  the 
initial  conditions  given  by  Eq.  £0.  and  equals 

,(«■»>  '«s/If(|-n-n.|=  -  «'» 

Here.  F  (  )  is  the  hypergeometric  function.  If  N  is  a  positive  integer,  the 
hypergennietric  function  in  Eq.  2 1  reduces  to  a  polynomial  of  degree  N  -  I. 

In  particular,  for  N  a  positive  integer 


N-l  , 

F(1-N.N.  I;r)  -  I  ♦  £  - flT'ii) - .(-*)'  (22) 

i  i  r(N-i)lr(i4i)r 

where  T  (  )  is  the  gamma  function  (riN  t  i)  *  nkn).  rd)  =  l.  r(i/2)  =  •/»). 

Eqs.  18a  or  18b  can  then  be  integrated  to  give  r(u.P)  Thus,  for  N  an  Integer,  an 
analytic  solution  for  the  interaction  region  is  obtained.  N«»te,  fro/n  Eq.  19b, 
that  y  "  (2N  +  l)/2N  -  I)  so  that  Y  i.  «=-/  3,  7/5  ;  •  ■  for  N  =  1 , 2,  3,  •  •  • 

The  variation,  with  time,  of  the  fluid  properties  at  r  -  0  can  be  found 
from  Eq.  21.  At  r  =  0.  n  -  p  -  a  so  that  r.q.  21  relates  t  and  a  (or 
p.p.T)  for  t  —  I.  The  variation  of  p  with  t  ,  at  r  0,  is  plotted  in  Fig.  3 
for  y  -  5/3.  7/5.  9/7.  11/9.  The  asymptotic  solution,  for  long  times,  ia 
indicated  by  a  dashed  line  m  these  figures.  It  is  seen  that  the  flow  approaches 
the  asymptotic  solution  more  slowly  as  y  decreases. 
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Figure  3.  Variation  of  Density  with  Time,  at  r  r  0,  for  Plane  («  =  0) 
Expansion  of  Initially  Uniform  Cas  into  Vacuum. 


The  entire  flow  field,  auer  long  times,  can  also  be  described.  The 
velocity  field  is  given  by  v  -  r/t  (Eq.  4).  The  density  distribution,  and  other 
state  properties,  are  found  frym  Eq.  2  1.  Note  that  a  +  p  =  2a  approaches 
zero  as  t  —  oo  .  Hence  (1  -  a)  (1  -p)/2(a  +  p)— oo,  as  t-*oo,  provided  a  <  1 , 
and  only  the  last  term  need  be  retained  in  the  polynomial  expansion  for  F(  ) 
in  Eq.  22.  This  results  in  the  following  expression  for  p: 


Eq.  23  is  valid  for  N  an  integer,  t  --  oo  and  a  <  1. 

After  long  times,  the  characteristic  line  a  -  1  propagates  with  velocity 
2 / ( y -  1 )  (from  Eq.  16a  with  a  -  1,  a  -  0)  and  is  thus  parallel  to  the  leading 
edge  of  the  simple  wave  (p  -  -1  in  Fig.  2).  With  increasing  time,  the  mass 
contained  between  these  lines  becomes  negligible,  (i.e.  the  simple  wave  re¬ 
gion  becomes  negligible)  and  the  leading  edge  of  the  expansion  can  be  approx¬ 
imated  by  R  2 1  /  (  y  -  1 )  -  Again  defining  p  5  r/R,  Eq  .  23  then  becomes 


pt  = 


(y-1)  r (Nt  1/2) 

-spr  TfN) 


3-y 


(24) 


which  inay  be  considered  valid  for  t  -  oo  and  0  <L  p  A  1.  Eq.  24  is  quite 
similar  in  form  to  Eq.  1  3b  with  it  -0, 

The  asymptotic  solution  noted  in  Fig.  3  is  obtained  from  Eq.  24  with 

p  ~  0 

Approximate  Asymptotic  Solution  for  Uniform  Initial  Density. 

Eq.  24  gives  the  asymptotic  nature  of  the  flow,  as  t  —  oo,  for  a  plane 
expansion  (cr  0)  of  an  initially  uniform  gas  cloud.  It  would  be  of  interest 
to  develop  corresponding  expressions  for  c  ylindric  al  and  spherical  expansions 
of  an  initially  uniform  gas  cloud.  Approximate  expressions,  which  satisfy 
conservation  of  mass  and  energy,  will  now  be  developed. 
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Eqs.  1 3b  and  24  suggest  that  a  suitable  lor  m  for  p  ,  as  t  -*  oo  ,  is 


ptff+1  =  D(l-r!Z)B  (25) 

where  B  and  D  are  constants  and  r\  ~  r/R.  Since  an  initially  uniform  gas 
cloud  is  being  considered,  R  will  be  taken  as 

R  --—t  (26) 

V-l 

The  constants  B  and  D  can  be  found  from  i  onse rvation  of  mass  and  energy 
The  latter  can  be  expressed  as 


1 

(T+  1 


CT  I 

PR  cIr 


(27a) 


1 

y  (y  -  1 )  (<t+  1 ) 


v  cr  . 

T  R  dn 


(271>) 


Eq.  27a  equates  the  mass  at  t  0  to  the  mass  at  t  *  and  Eq  27b  equates 
the  internal  energy  at  t  -  0  to  the  kinetic  energy  at  t  oo.  (The  internal 
energy  -0  as  t  -*u>).  Substituting  Eqs  25  and  2b  in  Eqs.  27,  with  v  -  r/t. 
gives 


B  ( VJLi)  _  i  ( n  1 )  N  -  1 

V  "  1 


(28a) 


D 


,  ,  .T+l  EtBrl  +  ~) 

1  .  /Y-  1 1  2 


c H  1 


r  ( n  + 1 )  r  (ill) 


1 8b) 
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z 

For  cr  =  1  and  N  an  integer,  D  =  (y'*-l)/4.  Numerical  values  for  B  and  D 
are  given  in  Table  I  for  cr  =  0,  1,  2  and  y  =  5/3,  7/5,  9/7,  11/9  and  13/11, 

Eqs.  25  and  28  are  in  exact  agreement  with  Eq.  24  for  a  =  0  and  N  an 
integer,  Since  Eqs.  25  and  28  are  an  approximate  integral  method  solution, 
they  may  be  considered  valid  for  N  not  an  integer,  as  well  as  for  cr  =  0 ,  1,2. 
A  similar  solution  has  been  presented  by  Ref,  6,  among  a  number  of  other 
approximate  solutions. 

The  variation  of  density,  with  q,  is  given  in  Fig.  4  for  y  ~  7/5  and 
cr  =  0,  1.  2.  It  is  seen  that  as  cr  increases,  the  density  decreases  more 
rapidly  with  q.  The  variation  of  density  wdth  q  is  given  in  Fig.  5  for  cr  =  1 
and  y  =  5/3.  7/5  and  9/7.  The  density  decreases  more  rapidly  with  q  as  y 
decreases . 

Free  Molecular  Expansion 

If  the  mean  free  path^  of  the  molecules  in  a  gas  is  of  order  Ro,  or 
larger,  the  expansion  may  be  treated  as  a  free  molecule  flow.  That  is,  col¬ 
lisions  between  molecules  can  be  neglected  and  each  molecule  moves  with  a 
constant  velocity  throughout  the  expansion  This. case  has  been  discussed  in 
Ref.  4  for  various  initial  density  distributions.  The  free  molecular  expansion 
of  initially  uniform  gas  clouds  is  discussed  herein  for  t  -*•  oo  and  the  results 
are  compared  with  those  from  continuum  theory. 

Consider  an  elemental  mass  m  of  particles  to  be  located  at  r  at 

o  o 

t  -  0  and  to  have  a  Maxwellian  distribution  of  velocities.  The  mass  of  par¬ 
ticles  having  velocities  between  v  and  v  +  Av  is  denoted  by  Attt  and  equals 


A  m 

o 

m 

o 


,  ,h.3/2 

4tt  (-)  e 


A  v 


(29) 


An  approximate  expression  for  the  mean  free  path  is  f  =  v~  /a 
v  is  kinematic  viscosity.  Free  molecular  flow  theory  applies  when 
i/R  =  v  /a  R  >0(1). 

O  O  O  O  — 


where 
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where  h  =  y/(2a  ).  Assume  that  the  container  which  confines  m  is  removed 
1  o  o 

at  t  =  0.  At  a  later  time  t  the  particles  in  the  velocity  range  v  to  v  +  Av 

will  be  contained  in  the  volume  bounded  by  r  =  r  +  vt  and  r  =  r  +  (v  +  Av)t. 

-2-3  -  °  ° 

This  volume  has  the  magnitude  AV  =  4tt  v  t  Av.  The  density  of  the  particles 
p  ( r ,  f)  ;  Am^ / AV  is  then 


r-r  2 
°\ 

...  .ill  _ 

ptf,  t)  ' _ 2k 


Eq.  30  gives  the  density  at  (r,t)  due  to  an  initial  elemental  mass  m^  at 
(r  ,  0).  The  expansion  of  a  gas,  with  an  arbitrary  initial  distribution  of  den¬ 
sity,  can  be  found  by  the  linear  superposition  of  the  elementary  solution  given 
by  Eq.  30  which  corresponds  to  a  spherical  (<j  -  2)  expansion  of  an  elemental 
three  dimensional  mass  source.  A  similar  expression  can  be  found  for  an 
elemental  line  source  (cr  -  1)  or  plane  source  (cr  -  0)  located  at  (r^,  0).  The 
result,  valid  for  all  o'  is 


where  rh^  equals  the  initial  elemental  mass  lor  cr  =  2 .  the  elemental  mass  per 
unit  length  for  or  1  and  the  elemental  mass  per  unit  surface  area  for  rr  =  0 
(assuming  a  symmetrical  expansion) 


Now,  consider  the  expansion  of  an  initially  uniform  gas  of  radius 
and  <r  =  0  1,2  After  long  times,  the  flow  may  be  considered  as  having 

originated  from  an  elemental  source.  The  solution  for  t  -*  oo  is  then  given 


by  Eq.  31  which  can  be  put  in  non-dimensional  form  by  introducing  r  =  r /U  , 

v  =  v/ a  ,  t  -  ta  /R  ,  etc.  In  the  present  case  define  q  =  (y-l)r/2t  so  that 
o  o  o  r 

r)  =  1  corresponds  to  the  leading  edge  for  the  corresponding  continuum  expan¬ 
sion.  (There  is  no  leading  edge  for  the  free  molecule  expansion  since  r  -*  oo 
as  p  -*  0).  Eq.  31  then  becomes 


>)t 


<T+  1 


=  O'  e 


-  B 1  T| 


(32) 


whe  re 


B'  -  2y/(y-1  )Z 

D'  -  yjz y  / tt  ;  y/2  ;  ^  yj---  f°r  a  ~  *  >  2 


Note  that  B'  is  independent  of  a .  Numerical  values  for  B'and  D'are  tabulated 
in  Table  I.  Mots  of  pt°"^  versus  q  are  presented  in  Figs.  4  and  5.  Although 
the  form  of  Eq.  32  is  different  from  that  of  Eq.  25,  the  variation  of  pta  with 
q  is  similar,  for  the  two  equations,  except  near  the  leading  edge  of  the  expan¬ 
sion.  (See  Figs  4  and  5.)  For  q  -  0,  Elq.  24  becomes  pt*7  ^  =  D  and  Eq.  32 
becomes  pt0^*  -  D'  so  that  D  and  D'  determine  the  variation  of  p  with  t  at 
q  =  0.  From  Table  I  it  is  seen  that  for  each  rr  and  y,  D  <  D'  so  that  the 
density  at  the  axis  decays  more  rapidly,  with  time,  for  continuum  flow.  The 
difference  in  the  expansion  rates  becomes  more,  pronounced  as  y  approaches  1. 

The  later  result  can  be  more  readily  understood  by  considering  the 
distribution  of  mass  within  the  free  molecular  expansion.  Let  c|j(q)  be  the 
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fraction  of  the  total  mass  which  is  contained  between  the  ordinate  r|  and  the 
axis.  Then,  using  equation  32, 


4  (q)  = 


r 

f  p'-1 

Jo 


dr 


oo 


<r  , 

p  r  dr 


for  (cr  =  0) 


(33a) 


=  1 


(<r  =  1) 


(33b) 


=  erf  (^/B7 q) 


2^  -B'q2 

"7^ 


(<T 


2} 


(33c) 


where  erf  (  )  is  the  error  function.  The  effective  leading  edge  of  the  free 
molecular  expansion  might  be  taken  to  correspond  to  4  (q)  =  0.  9, 0.  99  or  0.  999 
etc.  The  values  of  B''r)  which  correspond  to  4  (q)  =  0.  9,  0.  99  and  0.  999  are 
given  in  Table  II.  Let  %  2  be  taken  as  a  characteristic  value  defining  the 

effective  leading  edge  of  the  free  molecular  expansion  and  let  the  corresponding 
value  of  r  be  denoted  by  R  Then  R  ,f/t  *  Z\/Yh{.  This  may  be  compared 

with  R/t  =  2 / (y -  l )  for  the  continuum  case,  provided  y^  1  .  Hence,  the 

continuum  expansion  spreads  more  rapidly,  particularly  for  y  near  1. 


A  more  precise  comparison  can  be  made  by  computing  4  (q)  f°r  the  con¬ 
tinuum  flow.  This  can  be  done  in  closed  form  for  &  =  1 .  In  particular,  for 
<r  =  1  and  continuum  flow,  Eqs.  25  and  28a  show  4(9)  =  1  r  (l-q2)^7+^7 
For  a  given  value  of  «  *  1  -  4(h)  and  u  -  ]  ,  the  corresponding  values  of  q  for 
the  free  molecular  and  continuum  flows  are  related  by 

M  /  \  o  n  t  \f-Tru.  /  JTy  \J  1  -c  *  )/(Y+  1 )_  p-or  e«j  ancj  y  ^  l  (  t^e 

latter  gives  (R^^p,  ^  /(R^ ff )Cont  =  )  v~  1 n{  y  which  is  essentially  the 

comparison  given  in  the  body  of  the  report. 
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Table  II.  Effective  Leading  Edge  of  Free 
Molecular  Expansion 


4»(Tl) 

n/B'  n 

cr  =  0 

cr  -  1 

m 

0.  9 

1.  16 

1  52 

1.  77 

0.  99 

1.  82 

2  15 

2.  40 

0.  999 

2.  33 

2  63 

2.  85 

HYPERSONIC  JET 


As  noted  in  Ref.  5,  the  lateral  expansion  of  a  steady  hypersonic  jet, 
bounded  by  a  vacuum,  can  be  reduced  to  an  equivalent  unsteady  expansion. 
This  equivalance  is  now  discussed. 

Consider  a  nozzle  to  exhaust  a  hypersonic  stream  into  a  vacuum.  Let 

u  be  the  velocity  in  the  streamwise  x  direction,  u  be  the  value  of  u  at  the 
7  o 

nozzle  exit,  M  -  u  /a  and  6  be  a  c  harat  te  ns»u  slope  of  the  streamlines, 
o  o  o 

If  the  dependent  variables  have  the  lorm 

v/u  -  0(6)  p/p  =  0(1)  u /  u  =  1  +  0(6^)  (34) 

o  o  o 

then  the  continuity  and  v  momentum  equations  reduce  to  Eqs.  1  (neglecting 
terms  of  order  6^  compared  with  1)  with  t  replaced  by  x/u  .  These  equa¬ 
tions  can  be  solved  for  v  and  p  independent  of  the  x  momentum  equation, 
Bernoulli's  equation  can  then  be  used  to  find  u  (which  equals  u  to  order  6*") 
The  characteristic  slope  may  be  taken  equal  to  that  of  the  leading  edge  of  the 
expansion  after  long  limes 


6 


Y-l 


(35) 


The  above  solution  is  self  consistent  when  !*/»  1  (provided  y  is  not  nearly 
equal  to  1)  and  hence  applies  to  hypersonic  jets.  Thus,  the  previous  unsteady 
solutions  for  <r  =  0,  1  apply  to  two  dimensional  (<r  =  0)  and  axisymmetric 
(«r  =  1 )  hypersonic  jets  The  variable  t  in  the  former  is  replaced  by  x/uo  in 
the  latter.  If  a  non-dimensional  x  =  x/Rq  is  introduced,  then  t  is  replaced  by 
x/Mq  in  the  non-dimensional  equations. 

NUMERICAL  SOLUTIONS 

The  equations  of  motion  for  an  initially  unitorm  hypersonic  axisymmetric 
jet  exhausting  into  a  vacuum  have  been  integrated  numerically  for  =  10,  15 
V  =  1.4  and  Mo  =  10,  y  =  1 .22.  The  computations  were  performed  by 
A.  B.  Troesch  and  R.  F.  Kramer  using  an  existing  digital  computer  program. 
The  density’  along  the  axis  of  the  jet  was  found  as  a  function  of  distance  and  is 
plotted  in  Fig.  6.  The  approximate  asymptotic  solution,  previously  presented 
(Eqs  28)  is  also  noted  in  Fig.  6.  It  was  not  possible  to  continue  the  machine 
solution  ir.til  the  asymptotic  form  pIx/M^^  =  D  was  attained. 

Attempts  to  analytically  extrapolate  the  machine  results  to  very  large 
distances,  for  direct  comparison  with  the  approximate  theory,  also  met  with 
some  difficulty,  the  answer  depending  on  the  extrapolation  scheme  used.  For 
the  M  >  10.  15;  y  =  1.4  case  however,  some  confidence  can  be  placed  in  the 

ext  r  •;  ,olateu  value  D  -  0 .22  which  can  be  compared  to  the  approximate  value 
I)  ).  24  given  in  Table  1  for  y  =  1.4,  «r  =  l.  As  expected  the  results  for 
M  '  10.15  and  y  =  1.4  were  essentially  independent  of  Mach  number  and  are 
plotted  as  a  single  curve  in  Fig.  6b. 

Refe'*«*n..'  4  gives  numerical  results  for  the  density  at  the  center  of  an 
/Initially  uniform  spherical  ga3  cloud  expanding  to  a  vacuum  for  y  =  5/3.  The 
results  are  presented  by  a  plot  of  p  versus  t  which  yields  D  =  0.  38.  The 
latter  value  can  be  compared  with  the  approximate  value  D  =  0.432  given  in 
Table  I  for  <r  =  2,  y  =  5/3.  However,  it  should  be  noted  that  the  initial  por¬ 
tion  of  the  continuum  curve  of  p  versus  t  in  Ref.  4,  is  obviously  incorrect 
since  the  numerical  solution  for  p  versus  t  approaches  the  asymptotic  solution 
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DENSITY  (p) 


a 


DISTANCE  IX/M0  ) 

Variation  of  Density  wiD)  Distant  e 
along  the  Axis  (r  '  (1)  of  Axisvm- 
nietnc  steady  Hypersonic  Jet 
(M  -  10;  y  -  1-  22)  expanding  into 

vacuum.  Continuum  Flow. 


1).  VanatF  on  of  Density  with  Dis¬ 
tance  long  the  Axis  (r  =  0)  of 
Axis  y  rrammetric  steady  Hypersonic 
•Jet  (  =10,15;  y  =  !•  4)  expanding 

into  cuum.  Cofitinuum  Flow. 


from  below  rather  than  from  above  (compare  with  Fig.  6).  This  casts  some 
doubt  on  the  correctness  of  the  resulting  value  D  =  0.  38. 

In  general,  it  appears  that  the  value  of  D  given  by  the  approximate 
solution  (Eq.  28)  tends  to  be  somewhat  too  large  for  g-  =  1 , 2  although  it  is 
exact  for  u  -  0  and  N  an  integer.  The  two  numerical  solutions  noted  herein 
indicated  that  D  is  too  large  by  about  lO^o.  Further  numerical  solutions  are 
required  to  better  define  the  accurai y  of  the  approximate  solution. 


REFERENCES 


Sedov,  L.  I.  Similarity  and  Dimensional  Methods  in  Mechanics, 
Academic  Press,  New  York.  1959 

Keller,  J.  B.  "Spherical,  Cylindrical,  and  One  Dimensional  Gas  Flows", 
Quart.  Appl.  Math.  14,  1956. 

Courant,  R.  and  Friedricks,  K.  O.  Supersonic  Flow  and  Shock  Waves, 
Interscience,  New  York  1948 


Molmud,  P.  "Expansion  of  a  Rarified  Gas  Cloud  into  a  Vacuum", 
Physics  of  Fluids,  Vol.  3,  No.  3,  May-June  I960. 

Greifinger,  C.  and  Cole,  J.  "One  Dimensional  Expansion  of  a  Finite 
Mass  of  Gas  into  Vacuum",  RAND,  Rep.  P2008,  June  6,  I960. 

Stanyukovich,  K.  P.  ,  Unsteady  Motion  of  Continuous  Media, 
Pergamon  Press,  New  York  1  960 


i 


